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Abstract 



X/^ , This work explores the tensor and combinatorial constructs underlying the linearised 

^ ■ higher-order variational equations LVE^ of a generic autonomous system along a particular 

solution 0. The main result of this paper is a compact yet explicit and computationally 

P^ ' amenable form for said variational systems and their monodromy matrices. Alternatively, 

the same methods are useful to retrieve, and sometimes simplify, systems satisfied by the 

coefficients of the Taylor expansion of a formal first integral for a given dynamical system. 

^ \ This is done in preparation for further results within Ziglin-Morales-Ramis theory, specifically 

^^ ■ those of a constructive nature. 
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2 Linearised Higher Variational Equations 

1 Motivation and first definitions 

1.1 Introduction 

Integrability, an informal word for reasonably simple solvability, is an important problem in 
Dynamical Systems. Its opposite phenomenon, and specifically low predictability with respect 
to time, is usually summarised under the term chaos. If the system is Hamiltonian, as are 
most problems in Mechanics, the "chaos vs solvability" disjunctive is doubly advantageous. On 
one hand, it is amenable to the techniques of Symplectic Geometry. On the other, theory and 
empirics yield a specific, thus observable integrability condition: the existence of a precise amount 
of conserved quantities. 

The introduction of the algebraic approach by Ziglin, Morales-Ruiz and Ramis produced 
hallmark contributions to the study of the integrability of Hamiltonian systems [6, 22, 23, 30], 
essentially couched on a study of the invariants of a given matrix group, associated to a linear 
system: the first-order variational equations introduced in 1.2. 

A second step forward was carried out by Morales-Ruiz, Ramis and Simo ([24]) in order to 
extend the preceding theoretical framework to the Galois groups of the (linearised) higher-order 
variational equations along a particular solution. 

The second step described in the previous paragraph is the driving force behind this paper. 
A constructive version of the Morales-Ramis-Simo theorem was already started by Aparicio- 
Monforte and Weil in [2] and tangentially tackled from another viewpoint in [5] (see Section 5) 
and the present work aims at expanding this effort by offering a closed-form expression for the 
linearised higher variationals. May the reader bear in mind that at no stage in the results of this 
paper from Section 2 onwards is the system required to be Hamiltonian. 

1.2 Dynamical systems and variational equations 

In accordance with results succinctly described in 1.3 and thereafter, we need to observe the 
following convention outside of Sections 2 and 3: dependent and independent variables for all 
dynamical systems will be allowed to be complex. Any open set T C Pj, is an admissible domain 
for the time variable, embedded into the Riemann sphere to include t = oo as a valid singularity. 
Consider an autonomous dynamical system: 

z = X{z), where X : [/ C C" ^ C". (DS) 

Assume X is holomorphic. Basic mathematical objects are defined analogously to their real- 
valued counterparts: conserved quantities and foliations of solution curves. 

Definition 1.1. (DS) given, assume 



a) A first integral of (DS) is a function F : V D U ^ C constant along every solution of 

d 
' dxi 



(DS). Equivalently, such that DxF = 0, where Dx ■= J27=i-^i'^- 



b) For every z G U, let ip (t, z) he the unique solution of (DS) such that ip (0, z) = z, 
defined on a maximal open set I{z). Function ip : I (z) x U ^ C" thus defined is C^ 
and whenever y = (p(T,x) for some r G I {z), translation I {y) = I {z) — r holds and 
ip {t,y) = ip{T + t,z), for every t £ I {y)- 'P> is called the flow of (DS). 

Clarifying preliminary comments are in order whenever a particular solution (t) is considered: 

a) trivially, the partial derivatives -i^'-P {t, cf)) of the flow are multilinear forms of increasing 
order (or, alternatively, multidimensional matrices, see e.g. [17]) and may also be charac- 
terised as the blocks appearing in the Taylor expression of the flow along cf), minus the 



factorial denominators: 

^{t,z)^ip{t,cl)) + ^-^ — ^^"'^^+2!^ — i^-'/'} +3!^ — {^-0} +•••> (1) 

bracket notation summarising multilinear forms. 

b) each of these partial derivatives ^r 9? {t, cf)), inverse factorial unaccounted for, may also be 
characterised as satisfying an echeloned set of differential systems, depending on the previ- 
ous k — 1 partial derivatives and customarily called variational equations or systems. 
They are explicitly called higher-order whenever A; > 2. 



c) the variational system corresponding to A: = 1 is linear: 

G Mat„ (K) (VE^) 



dX 

Yi = AiYi, Ai{t):= — 



z=(t>{t) 



its principal fundamental matrix being the linear part of the flow along (p, and K = C {4>) 
being the smallest differential field containing C (t) and the solution. 

d) For fc > 2, however, the system is not linear, yet a linearised version may be found. The 
aim of the present paper is to do so with explicit formulae. 

1.3 Morales-Ramis-Ziglin theory and extensions 

Heuristics of all non-integrability results within the Ziglin-Morales-Ramis-Simo theoretical frame- 
work are firmly rooted in the following principle, expected to affect a widespread class of systems: 

// we assume general system (DS) "integrable" in some reasonable sense, then for 
every particular solution </>(t) of (DS) the differential system satisfied by each of 
the partial derivatives of the flow at cf) (t) must be also integrable in an accordingly 
reasonable sense. 

Any attempt at ad-hoc formulations of this heuristic principle for (DS) has an asset and a 
drawback: 

• there is a valid integrability axiom for linear systems (and thus, for (VE^)): that the 
identity component of an algebraic group attached to them, named the differential Galois 
group [22, 28], be solvable; 

• still, in order to transform this principle into a true conjecture it is necessary to clarify a 
notion of "integrability " for (DS). 

The latter item is cleared in the Hamiltonian case by the Liouville- Arnold Theorem establishing 
a sufficient condition for a system to admit, at least locally, a new set of variables rendering it 
integrable by quadratures. Said condition is the hypothesis on H in the following: 

Theorem 1.2 (Morales-Ruiz, Ramis, 2001). Let Xh be an n- degree- of -freedom Hamiltonian 
system having n independent first integrals in pairwise involution, defined on a neighborhood of 
an integral curve (p. Then, the identity component of the Galois group of the variational equations 
of H along (j) is an abelian group. D 

See [23, Cor. 8] or [22, Th. 4.1] for a precise statement and a proof. 

Theorem 1.3 (Morales- Ruiz, Ramis, Simo, 2005, [24, Th. 5]). Let H be as in the previous 
theorem. Let G^ be the differential Galois group of the k-th variational equations, fe > 1, and 
G := liniGfc the formal differential Galois group (inverse limit of the groups) of Xh along (p. 
Then, the identity components of the Galois groups Gk and G are abelian. D 
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Theorem 1.3 makes use of variational equations defined as in Section 1.2 and the language 
of jets, after proving said non-linear equations equivalent, in practicality, to any consistent 
linearised completion. Efforts towards a constructive version of this main Theorem, as well as 
the line of study described in Section 5, are hampered by a lack of consensus on the explicit 
block structure of said linearised completion. The present work, summarised in its main result 
(Proposition 4.5) aims at contributing to fill in this gap. Hence, outcomes will be purely restricted 
to symbolic calculus and do not constitute new results in the theoretical framework summarised 
above. 

Notation 1.4 (Multi-indices and lexicographic order). Part of the conventions listed below were 
already introduced in [5]. 

1. The modulus i = |i| of a multi-index i = (ii,...,i„) € Z"' is the sum of its entries. 
Multi-index addition and subtraction are defined entrywise as usual. 

2. We use the standard lexicographic order: {ii, . . . ,in) < (ji, . . . ,jn) means ii = ji, . . . , ik-i = 
jk-i and ik < jk for some A; > 1. 

3. Let F : U C. C" ^ C be a complex analytic function over the open set U . We define the 
lexicographically sifted differential of F of order m as the row vector 

Cam rp \ 

dx\'...dx'n"^ 7 

where ii + . . . + in = rn and entries are ordered as per <iex on multi-indices. 

4. We define 

Cn -\- k — 1 \ 
_ , L Dn± '■= dn,l+ dn,2 + ■ ■ ■ + dn±- (2) 

It is easy to check that the set of possible /c-ples of integers in { 1 , . . . , n} , or alternatively 
the number of homogeneous monomials of total degree n in k variables, has dn^k elements. 
Quantity Dn,k will become useful in Section 4.1 when LVEj, is introduced. 

Notation 1.5 (Multi-index binomials and multinomials). Given integers ki,...,kn > 0, we 
define the usual multinomial coefficient as 



ki-\ h kn] (ki -\ h kn] {ki-\ \- kn)l 



y ki,... ,kn J \ k J kilk2l---kn\ 

For a multi-index k G Z>q, deffiie k! := /ci! • • • A;„!. For any two such k, j, we define 

/k\ _^ kilk2l---knl ^ fki\ fk2\ (kn^ 

\V) ' PlW- ■ ■ ■ PJ {h - PlV- {k2 - P2y- ■ ■ ■ (kn - PnV- \Pl/ \P2/ \Pn^ 

and the multi-index counterpart to (3), 

/ki + • • • + k^\ ^ (ki + --- + kn)! 
i ki,...,k^ ] • ki!k2!---k„! ■ 



(3) 



(4) 



(5) 



2 Symmetric products and powers of finite matrices 

2.1 Definition and properties 

Let K he a field and V a if -vector space. Let us first recount the definition and the requisite 
existence and uniqueness results for the symmetric power of V. See [13, 14, 18] for details. 

Definition 2.1. An r symmetric tensor power of V is a vector space S, together with 
a symmetric multilinear map ip : V"^ := Vx .T. xV ^^ S satisfying the following universal 
property.- for every vector space W and every symmetric multilinear map f : V^ ^ W there is 
a unique induced linear map fQ-.S^W such that the following diagram commutes: 

V xVx .r. xV^^W 



/© 



s- 



In other words, isomorphism Hom (S, VF) = S{V",W) holds between the vector space of linear 
maps S ^ W and the vector space of symmetric multilinear maps V^ -^ W . 

Lemma 2.2. For any two symmetric r^^ powers {Si,(pi) and (S'2,V'2) ofV, an isomorphism 
Tp : Si ^ S2 exists such that (f>2 = V' ° V'l • 

Proposition 2.3. Given any vector space V and any r E N, 
a) a symmetric r power {S, (p) exists for V. We denote: 

• ViQ ■ ■ ■ QVr := ^{vi, . . . ,Vr) for every Vi, . . . ,Vr £ V, 

• Sym'^y , VQ .■'..QV orQ^V in place of S, 

k 

• V® := V Q ■ ■ ■ Q V for any vector v G V , and 

• i;®P := 17^ Pi ... 17®P", for any set of vectors Wi, . . . , w„ and multi-index p G Z>o- 

Conventions Sym V = V and Syui V = K arise naturally, 
h) Furthermore, {viQ ■ ■ ■ QVr : vi.,- ■ ■ ,Vr £ V} is a system of generators of SymJ'V. 

c) For any vector space W and multilinear map f : V^ ^ W , the linear map f^ induced by 
the universal property is defined on the set of generators ofSyufV as 

fo {Vl Q ■ ■ ■ Q Vr) = f {Vi, ■■■ ,Vr) . 

d) If dim V = n < 00 then every basis {ei, . . . , e„} of V induces a basis for Sym*"!/ ; 

{(ei .^.1. ei) (62 T?. 62) • • • (e„ .»■.". e„) : r^ > 0, |r| = r} ; (6) 
hence, dimSym^y = dn^r- 

In particular, symmetric products of vectors operate exactly like products of homogeneous poly- 
nomials in n variables. 

Remark 2.4. Sym*" may also be defined explicitly in terms of the tensor power 0^", delegating 
observation of a universal property on the latter and then taking quotients Sym^'F = 0^" V/ ~ 
modulo the equivalence relation vi<Si- ■ -01?^ ~ 'U(t(i)0' • •0'i'o-(r) for every a £ &r, thus equating, 
via isomorphism, Syioa^V to the subspace ({ej^ e^j • • • ej^ : 1 < ii < i2 < • • • < ir ^ n}) of 
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Given any ET- vector space W and two linear maps f,g:V^W, define 

h:VxV^Sym^W, /i K, ^2) := ^ [/ K) 5 M + / (1^2) <7 (^i)] • (7) 

Immediately bilinear and symmetric, it is granted a unique linear /iq : Sym V -^ Sym W, 
obviously defined /i© {vi V2) := h (t*!, 1)2), by the universal property. Let us write f Qg '■= /i©. 
It is easy to check that f Q g = g Q f and, given linear maps fi,gi : W ^ Wi, 

ifio f) Q {910 g) = ifiQ gi) o {f Q g) . 
A similar construction applies to the symmetric product of m, > 3 linear maps fi'-V^>- W: 
A • • • /m : Sym^F ^ Sym^W-, ^^ ...©«„ k^ -1 ^ /i K(i)) • • • /m K(™)) • («) 

Let us generalise the above symmetric product into one involving any two linear maps 

/ : Sym^i V -^ Sym^^ W, g : Sym^'^ V -^ Sym*^ W, ji ,J2,k,i2>0- 

Assume V and W finite-dimensional, V having basis {ei, . . . , e„}. We will use notation in 
Proposition 2.3. Defining the bilinear map 

(p{ui,U2) := U1QU2, Ui GSym'^V, i = 1,2, (9) 

we look forward to building a new symmetric bilinear function h in terms of / and g generalising 
(7), and proving there is a unique linear /iq completing the diagram 

Sym^iy X Sym-'^F — ^ Sym'^+'^W 

.^-^" (10) 

Sym^i+^'2y 

We want /i to be a symmetric, bilinear map depending on / and g and yielding coefficient 1 for 
all-round repeated vectors as in (7). Symmetric, multilinear h : l/'<Ji+J2 _> Sym*^"''*^!^ is easier 
to define, generalising (7): for any ui, . . . , Uj-^^^j^ G V, 

where 

Sn,j2 ■= W e 6ii+j2 : a (1) < • • • < a (ji) and a {ji + 1) < ■ ■ ■ < a (ji + J2)} ■ (12) 
Define 

{ipi X ip2){ui,...,Ujj^+j^) = {ipi{ui,...,Uj^) ,(/72('"ii + l,---,Wii+i2)) > 

(fi being the universal map of Sym-^^y; we intend the diagram involving the Cartesian product 

"J/XJ1+J2 



Sym^iy X Sym^'^y — ^ Sym'i+*2^ 



(13) 



to commute. Let Uji, . . . , Ui^ ^. be ji + J2 vectors in V, each an element of base {ei, . . . , e„}. 
We have 

((^1 X 992) (uii , . . . , tii,^+,, ) = {ui, • • • Wi^.^ , Ui^^^^ • • • w,^.^+^,^ ) 

and we may also split each set of vectors into copies of separate basis vectors: 

with IpI = ji and |q| = J2, and define k = p + q. The expression of (11) in these basis elements 
is now an immediate consequence of basic combinatorics: 



H 



ei ^.i.,ei,, 



i^n. 



, 072, ■ • •, ^n 



1 



E 



^ h > |P|=ii 






P2. 



P. 



|/(e«P)0 5fe®'^-^ 



leaving no option for (13) to commute but 



^(e®P,e®'i) 



1 



E 



^ h ) |P|=ii 



Pi +9i\ P2 + q2 

, Pi [ P2 ^ 



Pn + Qn \ J. f _CTP\ ^ ,V ^0p+q-P 



Pn 



|/(e«P)0 5(e«> 



(14) 



On the other hand, the universal property on the total symmetric product (Syra^^^^^V, (f 
yields a unique Hq such that Hq o ip = h, 

YXJ1+J2 




ipiXip2 




«l+*2 



w 



(15) 



Sym^iy X Sym-'^y 



Symii+i2y 



The fact 99 o (c^^ x 932) = <^ is immediate. And fixing f (and h) the uniqueness of Hq follows 
from construction: any other h, rendering (10) commutative would require the commutativity 
of the perimeter of (15), hence h, = Hq. 

Hence all we need to do is express f Q g := Hq in terms of its action on base elements (6) to 
obtain a simple, explicit form. 

Notation 2.5. When dealing with matrix sets, we will use super-indices and subindices in the 
following manner. 

1. The space of {i,j) -matrices Mat^„(i^) can either be defined by its underlying set, i.e. 
all dm,i X dnj matrices having entries in K , or as the vector space of linear maps between 
symmetric powers Homx (Sym-'y; Sym*VF) whenever V = K"^ and W = K^. 



2. It is clear from the above that Mat^' (K) is the set of all scalars a E K and Mat„' (K) 
(resp. Mat^'*^ {I^)) is made up of all row (resp. column) vectors whose entries are indexed 
by dn^k lexicographically ordered k-tuples. 

2. Reference to K may be dropped and notation may be abridged if dimensions are repeated 
or trivial, e.g. Mat^-' := Mat^^„, Mat^^„ := Mat^*_„, Mat„ := Mat^, etcetera. 

Checking product defined below renders diagrams (10) and (15) commutative is now 
immediate. 
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Definition 2.6 (Symmetric product of finite matrices). Let A £ Mat^'{^ {K), B £ Mat^;-^^ {K), 
i.e. linear maps A : Sym-^^i^"- -^ Sym^^K™- and B : Sym-'^K" -^ Sym^^K™. 

Given any multi-index k = {ki, . . . ,kn) G Z>q and |k| = ki + ■ ■ ■ + kn = ji + J2, 



{A B) (e^ . . . e^) = -i- J] (^ (^ef • • • e^) (iJe^-^^ • • • e^-^" 



(16) 



notation abused by removing to reduce space within basis elements (6), binomials as in (4) 
and summation taking place for specific multi-indices p, namely those such that 

IpI = ji and K Pi K ki, i = 1, . . . ,n. 

Tlie following is a mere exercise in induction: 

Lemma 2.7. The product of Ai,. . . ,Ar, recursively defined by 

^1 • • • Ar := (^1 • • • Ar-l) Ar, 
where At G Mat^':^', i = 1, . . . ,r, is expressed in terms of multinomials by 

(A,0...QA.)eQ'-= ,,,/+,, J2 L ^.)6^'^®^"' */ |k|=Ji + ---+jV, (17) 

sums obviously taken for pi + • • • + p^ = k and |pj| = ji, for every i = 1, . . . ,r. D 
Remarks 2.8. 

1. For an equivalent "non-monic" formulation of (16) (i.e. one for which entry i^i need not 
have coefficient 1 in its formal expression) using multi-indices in both columns and rows, 
see e.g. [8, 9, 10, 11]. 

r 

2. Notation in Proposition 2.3 extends to matrices: Sym''^ := A®"^ := ^ • • • ^. 

3. In the case of a square (1, l)-matrix A £ Mat„ {K), powers ®'' according to (16) and (17) 
are obviously consistent with multiple product (8), hence equal to established definitions 
for group morphism Syn/ : GL„ (V) -^ GL,„ (Sym*' {V)) in multilinear algebra textbooks 
such as the expression in terms of the permanent of A (e.g. [13, Th. 9.2]), or fi^(D • • • (§)A 
in multiple references such as [2, 5, 7]. 

Example 2.9. Given matrices A £ Matg' {K) and B £ Matg' {K), we may write them as 



A^l Aei \ Ae2 )^[ "^'^ "^'^ ) , B = ( Bef^ I Bei e^ I Bef 

V I ^ \ ^2,1 02,2 ' V i I I ^ 



/ 6l,l 6i,2 ^1,3 \ 

&2,1 &2,2 ^2,3 

&3,1 ^3,2 ^3,3 

V 64,1 &4,2 &4,3 / 



and it is immediate to check that the (4, 3) (hence four-column, five-row) matrix product 



-.03 



AQB={ {AQB){ef') (A Q B) [ef' Q e2) {Aq B) [ei Q e^') {AQB){e^ 
is equal to 

/ oi, 161,1 ^(fli, 261, 1 + 2ai, 161,2) 3(201,261,2 + 01,161,3) 01,26 

a2,2bl.l+2a2,lbl.2+ai,2b2,l+2ai,ib2,2 2a2,2hl,2+a2,lfcl,3+2ai,2fe2,2+ai.i&2,3 



1,3 



„ 7, : „ h "2,2"l,lT^U2,l"l,21-"l,2"2,l-t-''"l,lf2,2 ':"2,2U1 ,2^"2,1U1,3 -|-^"1.2U2,2-|-"1,1U2,3 „ J, I „ h 

O2.l0l,l + Oi,i02,l ^ n ' O2, 2D1. 3 + 01,202 



3 3 

a2,2b2.l+2a2,lb2.2+ai,2&3,l+2ai,ib3,2 2a2,2h2,2+a2,lb2,3+2ai,2fe3,2+ai.lb3,3 



:,3 



„ U I „ (> "2,2"2,l-t-i"2,l"2,2-M'l,2f3,l1-.!"l,ll'3,2 ':"2.2U2,2^"2,1 "2,3 1-^"1.2U3,2^"1.1U3,3 „ U I „ h 

O2,l02,l + Ol, 103,1 g ■ g ■ 02,202,3+01,203,3 

„ U \ „ h a2,2&3,l+2a2,lb3,2+ai,2b4,l+2ai 164 2 2a2,2b3,2+a2,lft3,3+2ai 2fc4,2+ai,l&4,3 „ h \ „ U 

02,163,1+01,104,1 3 3 02,203,3+01,204,3 

02,164,1 5(02,264,1+202,164,2) ^(202,264,2+02,164,3) 02,264,3 



The following is straightforward to prove from either direct application of the universal 
property or the analogous techniques used in [8, 10], and will not be delved into here: 

Proposition 2.10. For any matrices A, B, C , whenever products make sense, the following 
properties hold: 

a) AqB = BQA. 

b) {A + B)QC = AQC + BQC. 

c) {AqB)QC = AQ{BQC). 

d) {aA) Q B = a{AQB) for every a e K. 

e) If A is square and invertible, then {A~^) ' = (^®^j 
/; AQB = if and only if A = or B = 0. 

g) If A is a square (1, l)-matrix, then Avi Av2 • • • Avm = A^"^Vi • • • Vm- 

h) If V is a column vector, then {A B) d®(p+<?) = (^Av®^) {Bv®'^), for every p,q G Z^q- 

i) If rank {A) = r then rank {A®"") = dr,m and det A®"^ = (det A)V « i . D 
The next two results are immediate as well: 

Lemma 2.11. For any two matrices A £ Mat^-' and B G Mat^'^ and vectors vi, . . . ,Vj^q G V, 
the following holds, Sj^q defined as in (12); 

{AQB){viQ---Q v,+,) = -^ J2 ^ (^-(1) ® • • • ® ^-0)) ® ^ Ko+i) • • • ^^-0+'?)) • (18) 

\ q I creSj,, 

Proof. Nothing but the universal property on (11) and diagram (15) with different notation. D 

Lemma 2.12. Let /c > 1 and {ei, . . . ,e„} be a basis of K^ . 

a) (see also [8, 10]) For any A e Matf^'^ and B e Mat^'*", and every vector v G Sym'^K"-, 

{vQA)B = vQ{AB). (19) 

b) If {ei, . . . , e„} is the canonical basis formed by columns in Id„, then 

n 

Y: (e„0ldf-i) (e^0ld«'^-i) =Idf+\ 

m=l 

Proof. a) It suffices to prove it for basis elements of Sym'^: for any k such that |k| = p, 

{v A) Be®^ = vQ {AB) e®^. (20) 

But this is immediate from equation (18) or the definition (16) of itself, 
b) Using the previous item and the associative property in Proposition 2.10, 

(era Idf -1) (el O Wf "') = e,n (e^ Id®'^-!) = (e^ el) Idf "^ 

it is immediate to check that 6^06^ is a square (2, 2)-matrix whose only non-zero element 
is 1 in position m,m-, hence 



n 



^ (e„ e 



m ^ •^m) 



Idf -1 = Idf Idf ^1 = Idf +1. 

D 
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2.2 More properties of 

We need to generalise some of the above properties in Proposition 2.10 for later purposes. 
Matrices will not be necessarily square unless specifically defined as such throughout this Section. 

Lemma 2.13. Given square matrices A,Bg Mat„' ' and matrices Xi G Mat„'-'% i = 1,2, 

{A B) (Xi X2) = ^ {AXi BX2 + BXi AX2) , (21) 

and in general for any square Ai, . . . , A^ S Mat^''^ (K) and any Xi E Mat^'-^' {K), i = 1, . . . ,m, 

Cm \ / m \ 1 m 

o^Mox^] = ^J:oM)X^■ (22) 

j=i / \j=i / '"• creek i=i 

Proof. All it takes for ?7i = 2 is applying the universal property on either (7) (replacing V with 
Syni K"^) or (14) (replacing ji and J2 by k) on the product of A i? by each of the columns 
of Xi X2: indeed, for every i such that |i| = ji + J2, Xie®^ and X2e®'~P are both vectors of 
Sym K"^ whenever |i| = j'l, hence 

{A B) [(Xie®P) (X2e®'-P)] =^[A (X^e'^P) B (XacO'^P) + B (Xie®?) A (X^e'^'-P)] , 

and attachment of fj^^+j^-. J2\p\=ji (p) to both sides of the equation, along with (16), yields (21). 

Equally true for arbitrary m, using the universal property on the multiple product (8), now 
expressed on vi,...,Vm as (Oj ^i) (Oj ^'^) = 7^ E<7g6™ Oi^i^(T(i), replacing each Vi by the 
corresponding product 0[^]^ ^je®P' and attaching (^^'^■+^.'") ^ ( ^ ) as in (17). D 

Lemma 2.14. Given A G Mat^'-' and Xi, . . . , Xm such that Xi G Mat^'"^', assuming 1 < j < m, 
[a Id®"-^) (O x] = J2 [a{x,,Q...QX,^)]q O ^s. (23) 






\i=l / l^ii<---<ij^m s^ii,...,i 



Proof Defining B := Id®"'"^' G Mat;f~J''"~^' and Vi := Xie®P« where |pi| = Qi, i = 1, . . . ,m, 
equation (18) becomes 

{AQld^"'-^)QXre^P'- = -- Y. ^(^ne®Pn0.--0X,^.e®P'.)0 Q ^.e®^^(24) 

Attach C^^^ ^Qm\ y-> I Wq \yQi\^ sides of the equation and let k be any ordered 

multi-index having modulus qi + ■ ■ ■ + qm- The left-hand side becomes (A Id®™"-') (O™ i Xi). 
Multiplying the right-hand side of (23) by e and expressing the result in terms of each of its 
symmetric product factors, the product of the resulting five binomial and multinomial coefficients 
is equal to {"^ ^'^'") and the span of the multi-indices is precisely that of those in the right- 
hand side of (24) once embedded into its sum. D 

An immediate consequence of either Lemma 2.13 or Lemma 2.14 is 

Corollary 2.15. Given a square matrix A G Mat„' and Xi, . . . ,Xm such that Xi G Mat^'-'*, 

(m \ 1 m 

OX^\=-J2iAXi)Q[x^Q---QXiQ---QXm). □ (25) 

4 = 1 / j = l 
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Lemma 2.16. Given square matrix X G Mat„' , any vector v £ K^ and r > 1, 

{Xv Id®^) X®^ = X®''+^ {v Id®'') . (26) 

Proof. The first step is proving (^Xv Id®*^) X®*" = Xv X®*". This is immediate from (19) 
but let us elaborate on the proof for the sake of clarity and illustration. The fact Id®'' = Id^^^^ 
simplifies some steps. Equation (18) iov i = 1, j = 0, p = q = r yields, for every set of vectors 
vi,...,Vr e K", 

{Xv Id®'') {viQ---QVr)=XvQ{viQ---QVr). (27) 

Set Vi := Xe®P* G K"' for any given pj, slight notation abuse notwithstanding since |pi| = 1. 
Then (27) becomes 

{Xv Id®^) (Xe®Pi • • • Xe®P'-) = Xv Q Xe®Pi • • • Xe®P^ (28) 

Hence 

{Xv Id®'') X®^e®'' = {Xv Id®'^) ^ E ( "^ ) (^^®^' ® ■ ■ ■ ® ^e®P'-) 

^- pi,...,p,. \Pl' • • • 'Pr/ 

= i y ( ^ Ww Xe®Pi • • • Xe®P'- 

^!p±P API" ••'?-/ 



XvQ^ y I ^ 1 Xe®Pi • • • Xe®P'- 

^!pitrpApi"-"P^7 



= XvQX®''e®^ (29) 

sum multi-indices pj adding up to k and having successive moduli ji as always. 
Equation (20) and application of Proposition 2.10 imply 

X®''+^ {v Id®") e®^ = X®''+i (v e®'^) = Xv Q X®^^ef^' • • • X®^^e®^\ 

precisely (29). D 

If {K,d) is a differential field [28], i.e. 9 : i^ ^ X is a derivation, a linear map satisfying 
the Leibniz rule d{ab) = bd {a) + ad{b), this extends entrywise to matrices, d{aij) := {daij) 
and the Leibniz rule applies to 0. This is immediate and well-known for square matrices, but 
we are in a more general case: 

Lemma 2.17. For any given X G Mat^^i'^^ {K) and Y € Mat^2,i2 (j;^)^ 

d{XQY)=d{X)QY + XQd{Y). (30) 

Proof. It suffices, from expression (16), to check it true for symmetric products of vectors u,v ^ 
Sym ^^ '2 {K"') but this is as trivial as for homogeneous polynomials in n arbitrary unknowns 
K[Ei,...,En] in virtue of Proposition 2.3 or isomorphism Sym''i+'=2 (y*) ^ ^^^i+fe (y^ ^) ^nd 
can be gleaned from any of references mentioned in Section 2. D 

Although the next result will be rendered academic by simplified expressions in Section 4.1, 
it is worth writing for the sake of clarifying certain routinely-appearing matrix blocks a bit 
further. 
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Lemma 2.18. 

a) If Y is a square n x n matrix having entries in K and Y = AY , then 

-Sym^y = k(AQ Sym^^-^ (Id„)) Sym'^y. 

h) IfXe Mat^'-'i (K) and Y G Mat;^'-'^ {K) satisfy systems 

X = AX + Bi, Y = AY + B2, Ae Mat^ (K) , B, G Mat^'^'^ (K) 

then the symmetric product of these matrices satisfies the linear system 

^ (X y) = 2 (^ idd,^ ,) (X y) + (Bi y + ^2 X) . 

c) More generally, if Xi,Bi G Mat^'-'' {K), i = 1,... ,m, A e Mat^ and 

Xi = AXi + Bi, i = l,... ,m, 

then 



d_ 
It 



m m 



Ox, = m^40ldf™-MO^^ + E^^©O^r 



Proof. a) Immediate upon application of commutativity and (22) or (25) to 



(31) 



(32) 



(33) 



— (y0---0y =^y0---0 y q---qy = J2yq---q ay q---qy. 



dt 



b) Follows from Lemma 2.17 and the commutative and distributive properties of 0. 

c) We use induction, m = 2 being the previous item. Leibniz rule (30) and associativity imply 



d 



m—l 



/m— 1 



Ox,= [- C) xA Q x^+ [Q) xA Q x„ 



■ 1 .dt.^ , 

1=1 \ 2 = 1 / 

wherein induction hypothesis implies 

m— 1 m—\ 



.4 = 1 



i=l i=l 



jV* 



m— 1 



X^ + {AX^ + b^)qQ Xi, 



i=l 



which is equal in virtue of the distributive property and (25) to 



'm—l 



Li=l 



^ ^Xi (Xi • • • Xi • • • X^_i 
m—l m—l 

+AXm qQ Xi+B^Q Q Xi 

i=l i=l 

m 

= J2 ^^i Q(XiQ---QXiQ---QX„ 

.i=l 

further application of (25) ending the proof. 



m—l 



QX^+Y,BiQQ)XjQX„ 

4 = 1 j^/^i 



i=l jyti 



D 
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Remark 2.19. Albeit not explicitly as in (31), the matrix proven equal to k (AQ Id® 1 has 

appeared in numerous prior references (e.g. [2, 3, 4, 5, 7]) whenever an equation for Sym such 
as (31) arises, has been sometimes labelled sym^A and has been consistently called symmetric 
power of A in the sense of Lie algebras, its Lie group counterpart summarily standing therein 
for ®'^ as defined in this paper. 

3 Symmetric products and exponentials of infinite matrices 

The next step towards a compact form to linearised higher variationals is assembling the matrix 
blocks gleaned from Remark 2.19 together into a compact matrix whenever dealing with different 
blocks Yi,Y2, . . . satisfying different differential equations. 



3.1 Products and exponentials 

Of the myriad ways to note a set of infinite matrices, we may need one taking finite submatrix 
orders into account. Alternatively, of all the ways in which to write a if -algebra S, a need may 
arise to express it whenever possible S = Sym (V) := 0/^.>o Sym (V) for a given vector space. 

Notation 3.1. From now on Mat"'™ {K) denotes the set of all block matrices 

A = (^^j).,,;,o ' Aj ■■ Sym^if ™ ^ Sym^K\ 

hence Aij € M^^^xd^,, (K) = Mat j;f^ (if) ; 
/ •. 



A 



\ 



13,3 



i2,3 



Al,3 



A0.3 



A 



3,2 



A 



2,2 



Ai, 



A, 



0:2 



^3,1 



A2,l 



Ai. 



Ao,i 



A3.0 



A2.Q 

Auo 
Aofi ) 



We write Mat (if) := Mat"''^ (if) if the context allows for it, i.e. the value of n is unambiguous. 
Conversely, the set of matrices Mat^-'^ (if) is identified as a subset of Maf^'™" (if) by identi- 
fying every such block Aij with an element of Mat"''" (if) all of whose blocks are zero save for, 
perhaps, Aij. 

We define the following product on Mat"'™" (if). For a formulation yielding the same results 
see [11, p. 2]. 
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Definition 3.2. We define, for every given A,B £ Mat"'"" (K), matrix AqB = C G Mat"'™ {K) 
where for every given i,j, 



C, 



«j 



/ . ( • j^hji ® ^i~ii,j-ji- 



(34) 



^ il ^ i 
^ Jl ^ 3 



Same as always, ® will stand for powers built with this product. 

Example 3.3. Matrix AQ B takes the following form in its simplest echelons: 



/ 



ApfiBi^i + Aq^i Q Bi^o + Aifi Bq^i + Ai^iBp^o 



^o,o-Bo,i + ^0,1-^0,0 



^0,0^2,0 + ^1,0 Bifi + ^2, 0^0,0 



^0,0^0,0 



7 



and coefficients other than 1 will start appearing in further block rows and columns. We split 
row ^0 s-iid column g,* from the rest of the matrix for further clarity. 

The following is immediate and part of it has already been mentioned before, e.g. [10]: 

Lemma 3.4. Product on Mat"'*" is associative and commutative, and (Mat {K) , +, 0) is an 
integral domain as well as a K-algehra if endowed with the usual product by scalars in K . Its 
identity element is 1 made up of zero blocks except for o,o which is equal to !/<. D 

Definition 3.5. For every matrix A £ Mat {K) we define the formal power series 

expQ A:=l + A®' + ^A^' + ••• = £ ^A^\ 

1=0 

Whenever A has all but a finite distinguished submatrix Aj^k equal to zero (e.g. Examples 3.1 
below or Lemma 3.11), the abuse of notation expg Aj ^ = expg A will be customary. 

See also [10]. The fact is commutative saves us the trouble of having to check matrix 
commutation in the properties below, whose proof works exactly like that of scalar exponentials: 

Lemma 3.6. 

a) For every two A,B e Mat"-'"*, expQ {A + B) = expQ A expg B. 

b) For every Y £ Mat"''™' and any derivation 5 : K ^ K defined on field K , then 

5 expQ Y = {5Y) exp^ Y. 

c) ([8, Corollary 3]) For every two square matrices A,B£ Mat^' , expg AB = expg AexjpQ B. 

d) In particular, for every invertible square matrix A G Mat„' , expg A~^ = (expg A) . D 
Examples 3.7. 

1. Let A G Mat {K) such that all blocks are zero except for i^i: 

: \ 

Ma 



A 






^ 
^J 
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For 2^®^) expressions (34) containing Af-^ are those for which i = j = 2, hence 



l + A + -A®^ 
2 



/ ._ 




. \ 




^1,1 






Ai,i 




\ •• 




1 / 



and the pattern is clear in general: expg ^i^i = expg A = diag ( • • • , A®i , • • • , ^i,i, 1 

2. For row or column vectors this expression is even simpler. If the only non-zero block in A 
is a row vector, ^^^o = x = (^xi,..., Xd„j^j G Matn,fc,o (K), 

\ 



A 



/ . 








• ••• 





\ ■■ 


• X ■■■ 






the only expression (34) not automatically zero in {A A) is {A A)q 2k ~ i fc )^o,fc0^o,fc- 
Recursively, the only expression not automatically zero in A®^ is 



(a A^-^) 
For instance, if k = 1, 



^jk\ {j-l)k' 
. k ' \ k , 



2^ ] aGj 



(35) 



expQ 



A 






/ . 





























\ ■■ 


^®4 


^03 


x^'' 


X 


1 ^ 



3. This does not apply mutatis mutandis to matrices whose only non-trivial blocks are in the 
o,fc column. The only non-trivial block in A'^^ is jkfi whose expression is summarised in 
switching row and column indices and expunging binomials from (35). For fc = 1, we have 



expQ X = exp0 



/ . 




. \ 




• 







• 


X 


\ ■■ 


• 


; 



... 


J!^®' 


... 


^^®' 


... 


X 


... 


1 , 



A fourth example, namely matrices A G Mat (K) equal to zero save for block row i^^, fc > 1 
(see (37)), deserves special attention in the forthcoming Sections. Let us first fix conventions: 

Notation 3.8. For every set of indices satisfying 1 < ii < • • • < i^ and ii + ■ ■ ■ + ir = k, c^^ ^^ 
is defined as the amount of totally ordered partitions of a set of k elements among subsets of 
sizes ill ■ ■ ■ ) V- We will write c\ following i = (ii, . . . , v) cind omit super-index ^ if sum |i| is 
known beforehand. 



Remarks 3.9. 



1- (^ii,-,io = #h,Zk' following (61) below. Needless to say, En+...+i^=fcC: 

Stirling number of the second type ([1, §24.1.4]), and J2j=iJ2ii~{ hi =- 

/t*'^ Bell number [27, Vol 2, Ch. 3]. 



k 
k 



Bk, the 



{?}, the 
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2. Since each subset of size is is supposed to be ordered, we must divide the total amount by 
the orders of the corresponding symmetric groups, hence the explicit formula: 



^'^—'^ ni! •••«„! 



J112---JJ/ j {ii,...,ij) = {ki P-.K ki,k2 P?. k2,- ■ ■ ,km'^r^- km) , 



1 < fci < /C2 < • • • < k,n. 

Lemma 3.10. Let Y £ Mat {K) equal to zero outside of block row i^^, A: > 1; 



(36) 



Y :-- 



/ 








. \ 


















Y^ 


Y2 


Yi 





V 











/ 



K- G Mat 



1,1 



(37) 



Let Zr^s, s,r ^ 1, be the corresponding block in expg Y. Then, 

a) Row block r in expg Y is recursively obtained in terms of row blocks 1 and r — 1: 

1 *-'+i /A 

Zr,s = — / ^ \ ■ ]Yj Q Zj._ig_j. 

r p, \JJ 

In particular, Z^^r = Y®'^ and Z^^s = 0^^ ^^d^ s whenever r > s. 

b) For every m, r > 1 and any v € K^, 

(YiV Idf) Zr,r = Zr+l,r+l {v Id®'') . 

c) Using Notation 3.8 and (36), for every s ^ r 



(38) 



Zj'y^ 



E <,...,*.^n 0^*2 •••0^*.- 



(39) 



(40) 



^\^ Vtr 



d) Let A € Mat {K) defined with similar disposition as Y , its horizontal strip not necessarily 
at level i * ; 



A:-- 



/ : 






: \ 


... 











... ^3 


A2 


Ai 





... 











... 











V 








/ 



At G MafP'*. 



For every i, i > 1 and s ^ t + i, the following factorization holds: 

(t + i 



ll[M^tZt,)QZ,,s-j 



(^*0ldf) 



Z, 



t+i,s- 



(41) 



e) IfQ^ Mat" has only its square i^i block different from zero, then expg QY = (expg Q) (expg Y) . 
Proof a) Using (34) on ^ = y and S = y®^-\ 

^n.= ^ E (' MMi0^r-i,.-., =^ E (']Ys,Q{{r-l)\Zr-i,s-s,); 



hence, using the fact Zij = for i > j, (38) ensues. 
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b) Direct from (26) in Lemma 2.16. 

c) By induction. For s = 1, r can only be equal to 1 in order to have a non-zero block and 
Zi 1 = Yi = c\Yi. Assume (40) holds for all r smaller than or equal to s — 1. We have 



r ^ — ' \ 7 , 

j=l ^ ' iiH hir- l=s— j 



Summand redistribution renders the above equal to 

1 
r 



jl-\ \'jr=S 



(42) 



where, setting (ji, . . . , jr) = {ki .".i. ki, /c2 ."2. k2-,- ■ ■ ,km^^- km) , and defining kj := (ki ."■.\ ki 
and K.J :^ ( ki, . . . , Kj_i, fcj .'. . fcj, Kj_|_i, . . . , k^ I , 



Each of the summands in Cj^^...^,. is equal to 



'^4-^+ : 4?^ +•••+,; 4-^ 



fcl/"'^' ' 1/^2/ ^' 



(g-fcj)! 

s! fci!"iA:2!"2...fci_i!"i-ifei!"ifci+i!"»+i---fc^!""i 

/cj! (s - /cj)! ni!n2! • • • (n^ - 1)! • • • n^! 



TlmC 



'JlvJr 



hence the coefficient of l^i • • • 0ijv ™ (^2) is equal to ^ (ni + • • • + n^) '^ji,...,jr ~ '^jivjV 
d) Let us express the left-hand side in (41) in terms of (40): 



E 



rriiA \-mt=j 



mt 







fciH \-ki=s-j 



distributivity yielding it equal to 



E • E E cL„...,n^,C'..,fcJ^*(^™i0---0^-*)]0^fci0---0^fc- (43) 

j=t V"'/ rni,...,mt ki,...,ki 

The above multi-sums are indexed, respectively, by sets Jt.j and Ji,s-j, where 

Jr,c '■= {(^1) • • • ) ^r) G Z*" : 1 < ni < ri2 < • • • < rif. and ni -|- • • • -|- n,. = c} . 

The set of all ordered concatenations of index vectors in Jtj and Ji,s~j as j varies from t 
to s — i equals the complete set Jt+i,s- Conversely, for every multi-index 

n = (ni,... ,nt,nt+i,... ,nt+i) G Jt+i,s, 

consider the set of pairs of multi-indices m = {mi, . . . , mt) and k = [ki, . . . ,ki) whose 
ordered concatenation is n: 

It,i (n) = {(m, k) : a (mi, . . . ,mt,ki, . . . ,ki) = n for some a G &i+t} 

the terms in (43) indexed by Zt,i (n) are summed up in 



E 



(m,k)GXt,i(n) 



m 



Cm Ck [At (F„, • • • Ym,)] (Yfci • • • nj • 
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Let us discriminate among terms in the above sum. For every (m, k) G It^i (n), split m 
and k into copies of different integers: 

k = /ii.^l. /Xi,--- ,/Xg.^'. /xJ 



This obviously implies (equalling Mi or Ki multiplicities to zero whenever necessary) 

n = (^1 ^i+/^i /ii, • • • , /x^ax{p,g} ^— {--^It/^— {-"> /i^ax{p,,}) , (45) 

and 

max{p,g} irLax{p,5} 

the amount of permutations of n in (45) leaving m and k in (44) invariant is equal to 

(Mi + Ki\ (M2+K2\ /M^ax{p,,} + i^max{p,g}\ 

1^ Ml )\ M2 )'"[ M^,,|p,,| )- 

Multiplication of this product by Cn yields (writing r = max {p,q}) 

( Afi I ' \ Ma )"'\ Mr ) ' U _ Sl 



(Ml + Ki)\ (Ma + K2)l ■■■{Mr + Kr)\ MiIivTiI^uiI^^i+^i • • • Mr\Kr\^ir\^-+^- ■ 



(46) 



Let us now return to sum (43). Using multiplicities as in (44), the summand corresponding 
to a given (m, k) E Xt^i (n) has its coefficient equal to 



s 



Cm Ck 



I'Ml^lH \-Mp^p\ /--ftTl/ilH \-KqfJLq\ 



\\m\) \Mini + --- + MpfipJ Mi!---Mp! Kil ■ ■ ■ Kq\ 

which simplifies into (46). 

Hence (pc^c^"-' times [At (Fmi • • • 5^mJ] Ifci • • • Ifc^ equals c^ ,^ times ah per- 
mutations of the factors leaving these products invariant. This allows us to apply Lemma 
2.14 to At and Yq^ Yqi, := Ym, Q ■ ■ ■ QYmt QYk, Q ■ ■ ■ Q Y^r. 



\cLct' [^t{Ym,Q ■ ■ ■ QYm,)]QYu,Q- ■ -QYk, = ['^M (^i Id^) y0m0lbk. (47) 



i + V 



The fact every summand in (43) fits the same profile as the left-hand side in (47) allows 
us to factor (*^ ) ( yl^ Id^j out of the whole sum, namely ^i+t,s. 



e) Replacing each factor Yi. by QYi in (40) and applying Lemma 2.10 we obtain expg QY = 
(Zr^k) where 

hence matrix expg Y appears multiplied by diag (• • • , Q®^, Q^^, 1) = expg Q. 

a 

Lemma 3.11. Let A and Y be as in Lemma 3.10. Then, 

{A expQ Y) expQ Y = {AQ expg Id„) expg Y. (48) 
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Proof. Upon observation of (34), B := AQ expg Id„ S Mat {K) is defined recursively by 



Bi=Ai, 



Bh 



V 



(o)^. 



Let 



$1 = n 



^^ 




k>l. 



k>2. 



(49) 



V -^i.fc 



be the matrix formed by the first k row and column blocks in expg Y. Block row r of i? is 

Br := ( {,\)Ak+i-rQU'^'~' I C.zDAk-rQUr-' I ••• I (,:i)^i0ldr-' I I ••• I 

the row comprised of the first k blocks in A is written A^ = (Ai^fci ^i,fc-i ■ ■ ■ , ^i,i) and the first 
block column in $fc is Z'' := {Z^^k, -^fc-i,fc . • • , ^i.fc) . 
For every s > 1, block i^s in ^expg Y is equal to 



^s — AgZ — 2_^ Aj^j,sj 
hence for every r = 1, . . . ,k block r,k in (Aexpg y) expg Y is equal to 



fc-r+l 

E 

i=l V / 
which can be rewritten as 



\Xj Zr^i k-i 



k—r+l k—r+l 



fc-r+l 

E 

i=i 



/_^AtZt,j Q Zr-l^k-j, 



>t=l 



E E , (At^t,,)0^r-l,fc-„ 



t=l i=t 



vJ 



its innermore sum ostensibly calling for Lemma 3.10 reenacted with p = 1, s = k and i = r — 1. 
(41) indeed yields the above equal to 



fc-r+l / 
t=l ^ 

precisely block combination BrZ 



t + r -1 
r 



l^){AtQldT-')Z-, 



t+r-l,fc) 



D 



3.2 Application to power series 

Since polynomials and power series split into homogeneous components, Example 3.7(3) implies: 
Lemma 3.12. 

a) Let F €z K [[x]], x = (xi, . . . , Xn), be a formal series. Then there exists a set of row blocks 
Mp^ G Mat^*„ {K), i > such that F admits the expression F [x) = Mp expg X, where 



Mp :-- 



/ 







. \ 



i^ ••• 





) 



G Mat^'" {K) , X :-- 



( ■■ ■ 


. \ 


... 





... 


X 


iv ••• 


) 
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l,i 



b) If F = Fi X ■ ■ ■ X Fm is a vector power series, adequate Mp G Mat^*^ {K) render 



I 



F {x) = Mp expQ X 



where Mp :- 



M^^ M^^ 



\ ■■■ 



\ 



M 



1,0 



e Mat" 



J 



Following Definition 3.5, we write F (x) = MpexpQ x if it poses no clarity issue. D 

From the above Lemma it follows that every formal power series can be expressed in the 
form Mi^expQo;, where abusing notation once again 



Mp = Jp + M^'° := 



mI;^ m]^^ 



















+ 











M 



1,0 







(50) 



In other words: Mp equals the sum of two matrices with easily computable 0-exponentials: one 
following Example 3.7 (3) (same as a;) and one following pattern (37). This fact, Lemma 3.6, the 
fact (Mat {K) , +, 0) is an integral domain and the universal property of on finite products 
yield the following two results; see [8, 10] for a proof. 

Lemma 3.13. Given power series F = (Fi, . . . , Fm) and G = (Gi, . . . , Gp) in n and m inde- 
terminates, respectively, 

a) Ifn = m, Mpc = Mp Q Mq. 

h) expg F {x) = (expg Mp) (expg x) . 

c) MGoF = MGexpQMF. 

d) expg {Mg expg Mp) = (expg Mg) (expg Mp) . U 

Corollary 3.14. Let F [x) = (Fi, . . . , Fp) (xi, . . . , Xn) he a vector power series, y = F (x) and 



X 
Y 






R.,,x G Mat^'" (K) , Sy,Y G Mat^'^ {K) , 

be independent and dependent variable changes, which we assume admit formal inverse changes 

Rx,x e Mat"'^ (K) , SY,y G MafP'-P (K) . 

Then, the expression of F in the new variables, written in that in those old, is 

where y = F (x) = Mp^x,y^^P(T) x. D (51) 



X 

y 



Rx,xexpQX, 
SY,yexpQY, 



Mp^x,Y = Syy (expg Mp^^^y) expg Rx,x 



As was hinted at in [10, p. 5], the above result shows interesting light on the way finite- 
level transformations translate into transformations on Mat"'™". For a linear transformation of 
the independent variables x = BX, however, basic properties of expg are as useful as (51) in 
proving F admits the following expression in the new variable X (mind the effect of the first 
matrix, equal to zero save for block i.i which is equal to Id„, on the second one): 



F {X) = Id„ (expg Mp) (expe B) X = (jp + M"/ (exp^ B) expg X 



(52) 



This will be applied to first integrals of dynamical systems in Section 5. 
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4 Higher-order variational equations 

4.1 Structure 

Let us step back to what was said in Section 1.2. For each particular integral curve = 
{(p{t) : t G T C Pj,} of a given complex autonomous dynamical system (DS), the variational 
system VEj, for (DS) along cf) is satisfied by the partial derivatives ^r'/' {t, cf) {£)). 

Case k = 1 being trivial as shown in (VE,/,), the situation of interest is A; > 1. We will eschew 
formulations such as those in [24, eq (14)] in favour of a two-fold explicit expression: plain, sum- 
related expansion (63), and equally plain formulae (53), (64), (LVE^) and (VE^) using Linear 
Algebra to express multilinear maps. 

Notation 4.1. Define 4>{t) be a particular solution of (DS), K := C(0), Ai := X^'^' {(f)) and 
Yi := lex ( -^ifit, 4>) ) , and following Lemma 3. 10, let 

Zk-l,k 



$1 



Y, 



<^k 



\ -^i,fc 




k>2. 



(53) 



be formed by the first k block rows and columns in ^ = expg Y . Define A,Y (^ Mat (K) as in 
Lemma 3.10 with the above terms Ai, Yi as blocks. 

We also denote the canonical basis on K^ (meaning the set of columns o/Id„j by {ei, . . . , e„}. 



Lemma 4.2. In the hypotheses described in Notation 4-1, let k ^ 1. Then, 



dYk. 



>--i:1lf( 



m=l 



^m. ^ ^^n 



Qk-1 



(54) 



and for every m = 1, . . . ,n 
d 



Yk 



d 






0fc 



>^fc+i e„0ld^ 



■^r,fc+l ( ^m Id„ 



yie„0idf-Mz,_i,, 



d 



dZr. 



■Au 



Ak+i[YiemQldl 



&k 



(55) 

for every r ^ k, (56) 

(57) 



Proof. (54) is an immediate consequence of Lemma 2.12 and equation (55) which we are now 
going to prove. We have, for every given ordered multi-index i = (ii, . . . , if.), 



dYk 



-ei^ • • • e 



d 



aV 



Yk+iem Q ei-^ Q ■ ■ ■ Q Bi^ 



OZfyi CIZfYi OZi-^uZi2 ' ' ' CfZi^ 

The right-hand side in (55) is equal to this expression, too, by simple application of the same 
principle as in (20) in order to obtain ( e^ Id®^j e®' = e^ e®\ The effect of ^ on Aj is 
clear as well: chain rule implies 



dA, 



'-ei^Q---Qe. 



d 



dX 



E 



d^+^X 



dipr 






e^ e 



Di\ difr 



dzm^'^^ ""''''= dzm dzij^dzi2 ■ ■ ■ dzi^ ^^dzi^dzi^-'-dzi^dzrdzm ^^--«-^-^ V"' ""^ ^ dzm' 
which is equal, again using (20) in order to obtain gj^e,, e®' = (^7^6^ Id®**] e®', to 



^ j^ ^ / /^ 

^fc+i E {^r e®') 1^ = ^,+1 Y. i^^r e® 



r=l 



dZr, 



^^lV9^„ 



^ r=l 



n / n 



\dz„ 



Br Id^' 



0r \ „0i 
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hence to 



Au+i (^ Idf ) e«' = A,+i {Y,e^ Idf ) e«\ 



Let us prove (57) by induction over k. Assume the equation holds for all values smaller than 
k. Derivation of (38) and use of (55) and Leibniz rule (30) yields 



d 

dZrr 



■Zr-k — Sl + S2 — 5*3, 



where 



Si 



So 



I'-J^'fk' 



- E J y,+i e„0ld«^- 



]=1 V/ 



« ■ ;'s'e)'-- 



Zr-l,k-jj 

^r-i,fc-,+i(em0ldf-^) 

YiSra Idf -2) Zr^2,k-j 



Completion of ^3 with j = k — r + 2 and application of (41) with i = 1, t = r — 2, s = k, 
At = YiBm Id®''"^ and p = r -1 yields 

^3 = ^-^ {Yie,n Idf-i) Zr-i,k --(, ^ , }lYk-r+2 & [(^e^ Idf -2) ^,._2^^_2] . (58) 
r r \k — r + 2J 

The second term in the above expression can be written in the same manner as summands in 
5*2; indeed, using (26) and the fact Zr~2,r~2 = Y®^~ , 

(YiSm Idf ~2) Zr-2,r-2 = ^r-l,r-l (e„^ Idf-^) , 

hence the second term in (58) is 

namely the additional summand for j = A; — r + 2 in ^2 . 

A missing term ^ (YiSm 0ld®''~"^] ^r-i,fc needs to be accounted for in the first summand 
of (58). But it equals term j = for Si. Hence, index shift puts Si, S2 and the two extra terms 
from (58) together in a single sum: 



where 



^ k ' 



dz. 



Y, ( e„ Id®^'- 



J I ^m 



i=i 



QZr 



+ ( .15^,0 



Zr 



;-i,fc-j+i(em0ldf-^')]. (59) 



Let us check Sj is equal to ( i" ) (Yj ^r-i.fc+i-j) i^mQ Id® j . The columns of the latter matrix 
are of the following form, defining 1^ = 0, . . . , 1 , . . . , G Z" and whenever |k| = k: 



'k + 1' 



(y, Z,_i,,+i_,) 6®'^+!- = E ( ^ ^- ^" ) Y,e^' ^.-i,.-,+ie 

|i|=i V ' / 



0k+lm-i 
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and can be split into two sums depending on whether i = (ii, . . . , in) above satisfies im > 0: 






1 e^ e 



Ok-q 



|p|=i- 
precisely SjS^ as in (59). Hence 



dZ, 



r,k 



I'^^Vfc + i^ 



dz 



/r> / -^ 



K- 



i=i 



^r-i,fc+i-, (e,r^ Id®'^) - (yie„ Idf '-i) Z,_i,fc 



Zr,k+i (em Id®'^) - (Y^em Idf^i) Z,_i,fc. 



D 



Proposition 4.3 (First expHcit version of non-linearised VEj,). In the above hypotheses, 

y = ylexp0y; (VE<^) 

in other words, for every k ^ 1, 
d 



dt 



k k 

Yk = J2 ^jZk,j = T.^J E 4,...,i.^n QYi,Q---QYi^ 

j = l j = l ii-\ ^ij=k 



(VE^) 



Proof. Assume the result is true for k — 1, and let us prove it for k. That is, assume VE , 



k-l 



^k-1 



can be expressed in the form -^Yk-i = J2j=i ^j^j,k-i- We recall the entries in Ifc-i are partial 

derivatives 

(54) yields 



derivatives of (p{t,z), hence ^ = ^ on every entry, Schwarz Lemma applies and derivation of 



dt " ^,dt dz„, [''mOian ) 2^ 

m=l m=l 






dzm dt 



'm ^-^ -^^n y ) 



induction hypothesis and Leibniz rule render J^^fc equal to 



m=l 



dZm 



k-l 

/ , ^pZp^k-1 
p=l 






m=l 



k-l 



X^dAp dZp^k-i 

1^ ■Q^^P,k-i + A 
p=i ™ 



^ dz„ 



e^ 



Idf -1) ; 



equations (56) and (57) imply this is equal to Si + S2 — Ss, where 

n k-l 
m=l p=l 

52 = EE^P^M(^-0ldf-i)(e^0ldf-i), 

m=l p=l 

n k-l 

53 = EE^p(^i^-0ld0P"')^p-U-i(em0ldf-' 

m=l p=l 

Sum swapping in X]m Sp implies 

fc— 1 n k-l 

S2 = Y. ^pZp,k E (e- Idf -^) {el Idf -1) = 5] ApZp^, 

p=l m=l p=l 



(60) 
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Simple index shift on p and (39) render Si — S3 equal to 

n n 

J2 Ak (ne™ Id,f -1) Z,_i,fe_i (el © Idf -1) = Afe ^ (Fi Zfe_i.,_i) (e„ Idf "i) (e^ Idf -1) 

m— 1 m— 1 

which is equal to 

Afc (Fi Zfc_i,fc_i) ^ (e„ Idf -1) (e^ Id®^-i) = Afc (Fi Zfe_i,fc_i) = AfcZfe,^, 



m,=l 



D 



the missing summand in (60). 

The following is but a reformulation of the above result: 

Corollary 4.4 (Second explicit version of non-linearised VEj,). Let (p{t,cj)) = {(pi, ... ,(!)„) 
denote the flow of (DS). Let k ^ 1 be the order of the variational system. Given integers 
Ni, . . . , Nk > 1, r = 1, . . . , A; and mi,m2, • • • , ?i^^ ^ such that X]^=i "^-j = ^7 define: 

a) S]\fi,...,Nk '■= W (-^ij • • • ) ^k) '■ cr £ &k} and the set of partitions of {Ni, . . . , Nk} in ordered 
subsets of sizes mi, . . . , m^.: 

InI::S: ■= {(Ki, . . . , K,) e Sn„...,N, ■■ K, = {K,s, . . . , K„n,) , X.,1 < • • • < K^^r,^,} ; (61) 

b) and, using abridged notation J2ji,...jr ^'^ denote Z]ji=i Sj2=i ' ' ' Sjr=i' 






'Nu...,N. 



E E 



"Jr 



(Ki,...,KOe/ 



jTi]^ ,. .. ,mr 






dz 



Kr 



(62) 



Then, the order-k variational equation along 4> = {0(t)} is summarised in the following: 

k 



d 



E E ^;:;^; 



.,//tr 

,Affc ' 



i,A^i,...,iVfcG{l,...,n}, (63) 



dt dzN^dzN^ ■ ■ ■ dzN^ ^^^ mi,...,m, 

indices in X^mi ... m^ '^5'^^^ constrained by K mi < ?7i2 ^^ • • • ^^ Jti^' a?^c? m-i + • • • + m^ = /c. D 

In the previous Lemma we effectively settled the entries for lower n rows in A^y^k and the 

first n columns in $fc by virtue of (VE^). Let us now prove the result true for the rest of the 
matrices. 

Proposition 4.5 (Explicit version of LVEj). Still following Notation 4-1, the infinite system 



X = A 



LVE, 



X, 



Au/E^ ■.= Aq expg Id„, (LVE^) 

has ^ := expg Y as a solution matrix. Hence, for every fc > 1, 

a) the lower-triangular recursive Dn^k x Dn^k form for LVE^ is Y = Ary-p^kY, its system 
matrix being obtained from the first k row and column blocks of AlvEa •' 



^LVES 



f (,^)Ai0ld^ 

(fc'2)^2 0ld 



(o)^^ 



©fc-l 



ek-2 




(64) 
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b) and the principal fundamental matrix for LVE^ is ^k from ^ = expg Y in Notation 4-1- 
Proof. (48) in 3.11, (VE^) in Proposition 4.3, and item (b) in Lemma 3.6 imply 



exp0 Y = Y Q expg Y = (A expg Y) expg Y = (AQ expg Id„) expg Y. 
The rest follows from Lemma 3.10. 
Example 4.6. For instance, for k = 5 we have 

\ 

10^2 0ld®^ 4^i0ld®3 

^LVE5 = 

V 



D 



/ 5^1 Old®^ 

10^2 0ld®^ 4^i0ld®3 

lO^3 0ld®2 6^2 0ld®2 3^10 Id®2 

5^4 0ld„ 4^3 0ld„ 3A2 0ld„ 2^i Id„ 

^5 Ai A3 A2 Ai J 



and, using any of the equivalent expressions (38), (40), the fundamental matrix having Id/)^ 5 
as an initial condition is 



/ 



$, 



y05 

^03 



\ 



loyf •^ Y2 

lOyf ^ Ya + 1511 Y^^ 



y04 

6yf2 Y2 



Y 



03 



10^2 515 + 5^10^4 4^10^3 + 3^2 01^2 3^10^2 Y^^ 

V I5 I4 I3 I2 Yi J 



hence (VE . ) for /c = 5 can be expressed as 



Y5 = A1Y5+A2 {10Y2 QY3 + 5Yi y4)+^3 (ioyf2 QY3 + i5Yi Y^^) +A4 (loyf ^ ^2) +A5Y^ 

4.2 Explicit solution and monodromy matrices for LVE^ 

Let T C P^ be the domain for time variable t in (DS) and 7 C T a closed path based at point 
to £ T. Analytic continuation extends to polynomial functions, hence to symmetric products 
as seen in (16). Assume A; = 1. If Yi is a fundamental matrix of first-order (VEg,), analytic 
continuation along 7 yields 

Yi (to) ^^ Yi (to) • Mi,^, 

cont 

Mi^^ being the monodromy matrix ([31]) of (VE(j). Assume Yi := $1 is the principal fundamental 
matrix for (VE^), any other fundamental matrix ^1 recovered from ^1 = Yi^i (tg). 

Having computed Yi, the non-linearised second-order equation, after Proposition 4.3, is 

Y2 = A1Y2 + A2 • Sym2 (n) . (VE2) 

Following Proposition 4.5, linearised completion LVE^ has principal fundamental matrix 

$2 = 



05 



y02 



. ^2 Yi J- 
A particular solution Y2 of (VEj) is found via usual variation of constants: 

Y2 = YijY,~^A2Sym^{Yi), 
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which becomes a contour integral whenever time is taken along path 7: 

Y2 ^^ gi,2,7 := ^1,7 / n"'^2Sym2 ^y^^ ^ 

cont J7 

hence 

\ ^ , J M®2 



(65) 



Id„^..(.o)^^^"(^°^^- 







*^"°HQi; M.. 



cont \^ Yi (to) (3l,2,7 ^1 (to) Afi,^ y 

and [7] H^ Mi^^ is a group morphism vri (T, to) — > GL/)^^ - (C), hence for any fundamental matrix 



^2 (to) ^ ^2 (to) f ^^'^ ° 

cont \ '=t;l,2,7 -''^1,7 



therefore the monodromy of LVEj, along 7 will be 



M- 



02 



^2,7 



Qi,2,7 Mi,^ ; \^ Afi,^/^yf^^2n®' ^hn 



(66) 



Assume fc = 3. The principal fundamental matrix of LVE^ will be 



$, 



Sym3(yi) 
SFi Y2 Sym2 (Yi 



I3 



Y2 Yi 



and any other fundamental matrix can be expressed in the form ^3 = ^^C as usual. Let us now 
find a solution to 

Fs = AiY^ + ZA2Y1 QY2 + AsSym^ {Yi) , (67) 

Same as before, a particular solution Y3 of (67) is I3 = YiVa where 

1/3 = n^' (3^2^ I2 + ^sSym^ (Fi) 



(68) 



yielding a new contour integral if r G 7: 

Y3 ^^ Qi,3,7 := Mi,^ / yf 1 (3^2^! I2 + ^sSym^ (Y^)) dr. 

cont J-v ^ ' 



The remaining term of our monodromy matrix is a direct consequence of analytic continuation 
as performed on 3Yi Y2: 



= 3Yi (to) Y2 (to) -^ 3Mi ^ Qi,2,7 = 3Mi,^ f Mi,^ / Ff ^^2^®' 

cont \ J-y 

Our monodromy matrix is 



M 



3,7 



3Mi,^ Qi,2,7 M®2 



Q 



1,3,7 



Ql,2,7 Ml 



7 ^"1,7 



SMi^j Qi,2,7 
Ql,3,7 



M; 



2,7 



(69) 



The pattern is clear now. Assume we have computed solutions Yi, . . . , 1a;_i and performed 
continuation up to A; — 1: 



$ 



fc-i 



$fc_iMfe_i ^ := $fc_i 



/ <3A:-l,fc-l,7 

Qk-2,k-l,^ Qfc-2,fc-2,7 

Q2,fc-1,7 Q2,k-2,-y 

\ Ql,k-l,y Ql,k-2,'Y 



\ 



^2,2,7 

Ql,2,7 <5l,l,7 J 
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where 

Qr,s,'r ■= Yl c|^,...,j,.Qi,n,70Qi,i2,70 ••• 0<3i,v,7> s > r > 2. (70) 

Then, the fundamental matrix for (LVEi j will be expressed in the form (53), its lower left 
block Yfc being computable in terms of the blocks Z2^k, • • • > Zk,k above it (all of which involve 
Yi, . . . , y^-i) in virtue of (VE|): Yk = YiVk, which is continued into (5i,fc,7 := -^1,7 Xv ^k, where 

Vfc = Y"]" ^7=2 ^j^j,k- Upper terms Z2^ki • • • > -^fc.fc are continued into Q2,k, • • • > Qk,k as in (70), 
s replaced by k. It is clear we have proven the following: 

Lemma 4.7. The monodromy matrix ^^ o/LVE. along closed path 7 is composed by the first 
k row and column blocks in 

f ■■■ \ 

expg Q^ := expg ••• Qi,2,-y Qi,i,7 1 (71) 

y~^^ \o~) 

where (3i,i,7 := Mi^^, blocks above the bottom row are computed according to (70) and 

Qi,s,-y ■■= Mi,^ f Yf^ Yl ^i^i «' 2 < s < fc. □ (72) 

h j=2 

Hence it is clear the computation of a monodromy matrix follows a block order such as the 
one below, blocks in the bottom row requiring quadratures: 



7 








8 


4 




9 


5 


2 




10 


6 


3 


1 



(73) 



Computing the monodromy matrix is concomitant to computing the fundamental matrix, i.e. 
said bottom-row quadratures must be both indefinite (yielding terms Zi^g to be used in the 
computation of Zj^g in (72)) and contour integrals. See Section 6 for an example. 

We assume there are two generators [7] , [7] G vri (T; to); yielding two different matrices: 



7 



Q, 



7 



Q~. 



;7' ' ~*7 

Commutativity of monodromy matrices now admits simple, compact formulation: 

Proposition 4.8. Two monodromy matrices Affc_^ and Mj.^^ for LVEj, commute if, and only if, 
their previous blocks Mfc_i ,y,M^._i ;y commute and the additional properties hold 

k k 

/ J ^rjj-y^jjkj-y — / ^ ^r,j,y'^j,k,'jj jOr CVCry r — !,...,/£ i, 

j=r j=r 

matrices defined as in (70) and (72). D 
Remarks 4.9. 

a) The monodromy group of a linear system is contained in its differential Galois group 
(e.g. [28]). The motivation for the above Lemma and Proposition is to capitalise on this 
fact. This may in turn be a step towards future constructive incarnations of the Morales- 
Ramis-Simo Theorem 1.3. The main obstacle implementing Proposition 4.8, symbolico- 
computational issues aside, is the incertitude on whether -/Vf^^^ and M^^^ belong to the 

Zariski identity component Gal (LVE^ j ; a sufficient condition for arbitrary order is ful- 
filment at order 1, Mi^^, Mi^^y G Gal (VE,^)", itself an open problem in general. 
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b) All disquisitions and results on the variational jet in [20, 21] are referred to the lower n-row 
strip for commutators of these monodromies. More specifically: 

• what is called jet therein is lower strip Y in principal fundamental matrix $ = expg Y 
for infinite system (LVE^), and we will use this terminology in the following Section; 

• morphism properties imply monodromy matrices along path commutators are equal 
to monodromy matrix commutators: Mj^^^-i^-i^^^^ = M^^^M^j^^Mfc, ^3^^,71; 

• hence, the "jet commutation" properties in [20, 21] amount to lower strip Q^ 11 
(that is Y after passage along 72" 7^^ 7271) equalling Id^. 

Although [20, 21] clearly benefit from the use of automatic differentiation techniques (see 
also [19]), it may be argued that expressions such as those in (LVE^) provide for a fuller 
control of the general structure of the whole variational complex when it comes to sym- 
bolic computations, as well as a further check aid for the aforementioned techniques. See 
Subsection 6.1 for an example. 

5 First integrals and higher-order variational equations 

Let F : U CI C" — > C" be a holomorphic function and (p : I C C ^ U. Firstly, the flow if (t, z) 
of X admits, at least formally, Taylor expansion (1) along (f) which is expressible as 

V9(t,0 + O = <^ + n^ + ^>^2^®' + --- = + ^</,expQt (74) 

where J^ is the jet for flow ip{t, •) along cf), displayed as Y in (37) and defined in Notation 4.1 
- that is, the matrix whose ©-exponential <I> is a solution matrix for (LVE^). 

Secondly, the Taylor expansion of F along cf) can be written, cfr. [5, Lemma 2] and Notation 
1.4, as 

F (y + 0) = F (0) + ^ — (f(-) (0) , Sym"y ) . (75) 

m=l 

Basic scrutiny of Example 3.7(3), Lemma 3.12 and (50) trivially implies (75) can be expressed 

as F {y + (p) = Mp expg y, where 



/ ••• 

m| = J^ + F(o)(0):= ••• F(2)((/>) fW(0) 





fW(<^) I GMati'"(i^) 



\ ■■■ 

i.e. Jp is the jet or horizontal strip of lex-sifted partial derivatives of F at cj). 
Definition 5.1. We call 

Alyei ■.= -{AQ exp0 Id„) ^ , (LVEl) 



^ — Au\i-E*X, 



the adjoint or dual variational system of (DS) along (j). Same as in (LVE,^) and all throughout 
4-1, consideration of finite subsystems, namely the lowest Dn,k x Dn,k block, leads to specific 
notation (LVE^ j . 

The following is well-known for finite systems and immediate upon derivation of equation 
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Lemma 5.2. (^^ ) is a principal fundamental matrix of f LVEi 1 for every A; > 1. 
Hence, {^~^) , is a solution to (LVE^), where <I> = expg J^. D 

The following was proven in [24] and recounted in [5, Lemma 7], and may now be expressed 
in a simple, compact fashion: 

Lemma 5.3. Let F and cj) be a holomorphic first integral and a non-constant solution of (DS) 
respectively. Let V := Jp be the transposed jet of F along cf). Then, V is a solution of (LVE*). 

Proof. Let us recall formal expansion (74) and F (y) = JpexjpQy for every y G K^. Let tp = 
if{t,cf) + ^). We have, using Lemma 3.13, 



F {ip) = F {(f) + J0 expQ $) = M| expQ ( J^ exp^ $) = (^M| expQ J^j exp^ $, 
and F (ip) is supposed to be constant, hence applying (LVE^) and Lemma 3.13 



= i^Mf, expg J0 j expQ | = I Af | + M^^lve^ j expg J^ expg, | = I Af| + A/^^lve^ j expQ (■0-0) 



hence Mp + MpAu^^ = leading us to the final result after transposing both sides. 



D 



Compound the jet of field X, i.e. A in Notation 4.1 and Proposition 4.5, with a i_o term ^O) 
equal to X(o) = X (0) = 0: 



A:-- 



... 








... A2 


Ai 


Ao 


... 









^i:=X«(0)GMati'^(K). 



It is easy to check, via possibilities offered on ii and ji in (34), that the symmetric product of 
A with expQ Idn adds only a relatively minor addendum to ^lve,,, > namely a superdiagonal of 
blocks (*)^o Id®* G Mat^+^'*, i > 1, effectively rendering it block-Hessenberg: 

Alve^ := AQ expg Id„ = hm ^lve*= ' 

k 4> 

where, isolating ^lve'-' within ^lve* ^^ means of a solid line. 



A 



LVE 5 



Ao Id: 



0fc 



(5A,_.0idf 

QAk-lQidn 

Ak 


Ao0ldf-i 


Ao0ldf 




(^■-l)Afe'3 0ldf .. 

{''-') Ak-2Qldn ■■ 
Ak-1 




Ai Id„ 

^2 


Ao Id„ 


Ai 















Ao 
/ 



(76) 



(^)x("M'/')0idr 


\ 


(,^)xu)(0)0id--^ 


^LVE^i 

/ 


(,!jx(^M</')0idr"' 




(^-)x('=)(0)0id-'^ 



V 

Using the M^-Mk notation in [5], it is immediate to check that 



Mi 



Id®'^-! = Id®'^-' X (0) 



(77) 
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^LVE'= ~ -^fc-i ^o^ every /c > 1. A result in [5] using said notation is easier to prove 
in this setting. Indeed, the same reasoning underlying (55) applies to row vector F^^\ hence 
_g_^(fc) ^ ^(fc+i) J^g^ Q i^Qk\ ^ ^^^ following Lemma 2.12 



j,(k) ^ j,ik+i) J2 (e^ Idf ) 0„ = ^(^=+1) (0 Idf ) = f(^+i) [Ao Idf ) , 



m=l 



implying ( i^^*^^ j = (AqQ Id®*^) (i^C^+i)) ; placing all terms on one side, and observing Lemma 
5.3 and the transpose of expression (76), we obtain: 

Proposition 5.4 ([5, Th. 12]). Let F, cf), V be defined as in Lemma 5.3. Then A[y^ V = 0. D 

Hence, blocks in Vi, {V2, Vi) , (^3, V2, Vi) , . . . having all entries in the base field K and 
satisfying both equations in Proposition 5.4 and 5.3 are candidates for jet blocks F^^\ F^"^' , . . . 
of a formal first integral. These blocks belonging to the intersection of ker Aj ^. and the solution 

subspace Sol/^ (LVE^ I were called adm,issible solutions of the order-A; adjoint system in [5]. 

This takes us back to the end of Section 3.2. Consider gauge transformation ([2, 5, 6, 22]) 
X = PX transforming linear system ^ = ^1^ into equivalent 



E = P[A^]E: 

Using notation Yi = PXi, J^ = PX and item (e) in Lemma 3.10, we recover the result already 
seen in previous references, summarised in the extension of gauge transformations to higher 
dimensions via P®^: 

expg {X) = expQ {P^^J^) = expg P'^ expg J^ = diag (^- • • , (P^^)®^ , (P^^)®^ , P"\ Ij expg J^, 

and very simple application of properties seen so far extends the general structure of the gauge 
transformation to ^ = expg P~^ expg J^: 



^ = P 



A 



LVErf 



^ := expg P-^LVE^ expg P + { P~^ Q exp^ p-^ exp^ P ^. 



(78) 



The second term P expg P ^ J expg P inside the brackets can be simplified into the fol- 
lowing: 

,0fc-i 



diag . . . ,k 



p-'q{p~') 



P®^ . . . , 2 P-I P-M P®2, -P-^P, , (79) 



with P = — P ^PP ^. The above gauge transformation can be seen as the effect of transfor- 
mation z = PZ on the jet of (DS). Given a first integral F of the latter, we may always assume 



rl,0 



F {(f)) = 0, which implies Mp = and, as seen in (52) or in Lemma 3.6, 

Fp (Z) = Jf (expg P) expg Z. 
The jet of this formal series is 



Jpp = Jf (expg P) 





FW{(j))pe3 i7(2)(^)p02 i7{l)(0)p 















G Mat^'" (K) 
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and applying (78), Lemmae 5.3 and 5.4, and the identity 

I '■ — 

,0fc 



P 



-1 



P 



&k 



p-lfV^ 



we have just proven the following: 

Proposition 5.5. The transposed jet Vp := Jp in the new variables must satisfy 



Vp = -P 



A 



LVE. 



"" Vp, vIlve, (expg P-^S'vp = 0. 



D 



(80) 



The key importance in practical examples resides in the choice of the particular solution (p 
and the reduction matrix P, in order to render (80) easier (or more convenient) to solve than 
its unreduced counterparts. Lemma 5.3 and Proposition 5.4: see also [3, 4]. 



6 Example 

The dynamics of the Swinging Atwood Machine (SAM), summarised in the diagram below. 





are governed by Hamiltonian 



«4 



pj ^ {P2+RPI? 

Mt mqf 



+ gqi (M — m cos ^2) — gR{Mq2 — m sin (72) ; 



where Mt = M + m + 2/p/i?^, Ip is the pulley inertial momentum, qi = r, q2 = 9, pi = Pr, and 
P2 = Pe- We know ([26, Th. 7.5]) the following: 

Theorem 6.1. For every physically consistent value of the parameters, regardless of Ip and R, 
7i is meromorphically non-integrable. 

For SAM without massive pulleys, i.e. the limit case Ip = 0, R = and Mt = M + m: 



+ gqi (M — m cos 52) 





Ti-w = 


1 
" 2 


( Pi 

yM + m 


+ 


Pi 
mql 


the following holds: 












Theorem 6.2. 












1. ([15, Theorem 


1]) IfM> 


m and 







M , p(p+l) 

/i = — / ^p :^ 



vn ' p{p + 1) — 4 
for every p € Z,p > 2, then Hamiltonian Tiyj is non-integrable. 
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2. ([29, Equation (16)]) For p = 2, fj, = H2 = 3, X-h^, is integrahle and has the following first 
integral: 

T 2 ■ ( ■ QlQ2 \ , 2 2 2 2, Pl9lC - 2p2S 

I = (li<l2\qic i^sj+gq^sc = gqiC s + p2 — ^;^ , 

where c = cos ((72/2), s = sin ((72/2). 

3. ([20, Theorem 4]) The degenerate cases Hp, p > 2 in item 1 are non-integrahle. D 
Canonical transformation 



11 = Qi, Q2 = arccosQ2, Pi = Pi, P2 = --P2\/ 1 - Qh 
performed on Tiyu yields Hamiltonian 

2 \M + m mQi j 

Let us apply the techniques introduced in the previous Sections to the variational systems for 
this Hamiltonian. 



6.1 Monodromy matrices and integrability 

Consider the particular solution 



<A = (Qi,Q2,Pi,P2)= \- 



'{t-l)t g{m + M) {2t - 1) 5^"?, {t-l)t (Ci - 2Cit + t^) \ 

2 '"^' 2 ' 4 (Ci - 1) j 



(81) 



In the forthcoming calculations, any value of Ci different from or 1 will ensure the presence of 
logarithms in the fundamental matrix, and any value of Ci different from 1/2 will avoid division 
by zero. Choose, for instance, Ci = 1/3. Gauge transformation x = PX with 



/ 



P:-- 



\/M + m 



-1 










2(M+m)(l-3t)^ 

(M+m) (l+f)(3f--l) 
9{t~lff^ 



V 



gm(l+t) g(l5m(t-l)t(3f-l)'^+Af(l+t(16+15t(-9+i(37+27(-2+t)t))))) 

' 2(3t-l) 90(l-3t)^(t-l)t 







9grn{t-Vfj?_ 
2(l-3t)^ 



transforms (VE,j,) into the simplified system X = P [Ai] X where 



P[A 



I 



iG 



1 



9 \{t~if *■ 





%M{1-Zt) 
405m(t-l)4t4 



K 



{t-iY 



(82) 



A fundamental matrix for this system is 



* 



/ 1 






l+3t 

9(t-l)t 

1 





\ n 4M{§(2t-l)[l+5(t-l)t(6(t-l)t-l)]+t[3+5t(3-2t(ll+3t(2t-5)))]-20(t-l)^t^ In -^i} 

\ 405m(t-l)''t3 



t 



1 

l+3t 
'9(t-l)t 






1 / 



because P \A{\ commutes with its own indefinite integral - an interesting topic for further study 
is the possible relationship between said commutation and the level of reduction of the matrix 
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in the sense explained in [3], but this is not central to our present work. Our intention is to 
perform analytical continuation along two different paths, shown below: 




Along 7i, all terms outside of the diagonal in ^ vanish except for term (Mi^-y^)^ ^ which can be 



obtained from X = P \A-\\ X. by assuming X2 
and thus the contour integral 



1, X3 = 0, which yields 



8M(3f-l) 
405m(t-l)*t4 



+ X4 



(Ml, 



7lM,2 



71 



8M(3t - 1) _ 32iM7r 
mhm{t-\ft^ ~ 81m 



First-order monodromy M\^y^ along 72 is obtained analogously as Id4 plus sub-diagonal term 



(Ml,7. 



72 74,2 



72 



8M(3i - 1) 



32iMTT 
81m 



Second- and third- order monodromies are given by (65) and (68) respectively. The only change 
therein is replacing Ai, A2 and ^3 by their gauge transforms defined by (78) and (79) as the 
lower row block in 



A 



LVE? 



102 



^ 3 [gQ3 (Ai Idf ) +QQ> 
3QQ2(A2 0ld„)PQ3 



p03 



QQ2 (A^ id„) +QQQ 
p-lA2P®2 



p02 



\ 



^[^1] / 



where Q := P ^. Let Yi = ^. Following the order described in (73), we first compute upper 
indefinite block Y"® in the fundamental matrix, as well as the corresponding block M® in 
the monodromies. We then place indefinite block Y"® under the integral sign to obtain both 
the indefinite integral Yi JY{~ (P~^j42-P®^) Sym (Yi) defining Y2 and its contour counterpart 
defining <5i,2,7j- Using indefinite terms Yi and Y2 (resp. Mi^y. and (51,2,7,) we then build terms 
Yj®^ and 3Yi Y2 (resp. MfJ"^^ and 3Mi,-y^ Qi,2,7J and we finally plug Yf ^ and 3Yi Y2 into 
the corresponding contour integral for Qi,3,7i as in (68) (where ^3 is replaced by P~^A^P^^). In 
order to decide whether or not we also need indefinite integral Y3 for higher orders, let us check 
whether our third-order monodromy matrices commute. We have C := M^^^^^M^^^^ 
equal to zero except for the following terms: 



■^3,72^^3,71 



Q 



25600 
, 85293 



8960 
85293 



i) M\^ 



34,5 



g^m?{M + m) 



a 



/ 44400640 
I 209564901 



+ 



21766615859 2 
69854967 



34,11 



i M^vr' 



3^3 



C: 



/ 2220032 
I 2587221 



+ 



g'^m?{M + m) 

2775040 ■\ ^2 2 
7761663^ ' ^^^ " 



34,13 



g'^m?{M + m) 

This, coupled with the fact that all order-A; monodromy matrices built from Mi,^. belong to 
GalfLVE^j (since Mi,^^, being unipotent, belong to Gal(VE(^)° and all higher-order funda- 
mental matrices are obtained from quadratures), allows us to complement the non-integrability 
proof in [20] (see Theorem 6.2, item 3) using linearised variational equations. 
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Remark 6.3. Although we used monodromies in our argument, a similar reasoning applies to 
Galois groups by noting the presence of In (3t — 1) inY^ and Ini, In (t — 1) in Yi,Y2, Y3 and using 
the Picard-Vessiot extension C (i) (lni,ln (t — 1) ,ln (3t — 1)) | C (t) (lnt,ln (t — 1)) to glean the 
structure of a generic Galois group matrix on the fundamental matrix ^3 or its pre-gauge 
counterpart diag (P®'^,P®^,P) ^3; see e.g. [25]. 



6.2 Formal first integrals and admissible solutions 

Let us now try to apply gauge transforms to the adjoint system for the same Hamiltonian. We 
have a particular solution 



(f}={Ql,Q2,Pl,P2) 



-|(t-l)t,-l,-|(M + m)(2t-l) 



g^m 



{t-l)t' 



which corresponds to special case Ci = in (81), and gauge transformation x = PX with 
/ -7^T^ \ 



P:-- 



VM+m 



-VmT 

gm{l-t) 



m 



\'M+m 
g m(t-l)^ 

g{3m~M) _ gt/j^ _|_ 



m 











VM+m I 



12{t-l)y'M+m 

transforms (VE^) into the parameter-free, simplified system X = P [Ai] X where 



P[Ai 



( 






(t-i)^ {t-\Y 



{t-\f 



0\ 






The principal fundamental matrix for this system is 



/I 



^ 



{t-\Y 



1 





(i-t) 


1 








(i-2)f 1 
2(t-l)^ 



and P 



A 



LVES 



is computed as in the previous Subsection. In the following table, the first two 
columns correspond to order k and total solution space dimension -D4,fc. The latter two display 
the dimension of the subspace of rational solutions (i.e. those having their entries in the base 
field K = C (t)), and the subspace of those among the former satisfying (80), respectively: 



k 



dime SolP 



LVE5 



dime SoIa'-P 



LVE 



dime Soladm [P 



LVE 



1 


4 


4 


3 


2 


14 


14 


9 


3 


34 


32 


17 



Using solution (81) in Subsection 6.1 and gauge transformation P given in (82), however, 
dime Soli^P LVE . is considerably reduced and we obtain the following table, displaying 



smaller bounds on the amount of admissible solutions: 



dime SolP 



LVE5 



dime SoIkP 



LVE 



dimcSoladm (P 



LVE 



1 


4 


3 


2 


2 


14 


9 


5 


3 


34 


19 


9 
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